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We investigate thermodynamic and dynamical stability of a family of six-dimensional braneworld 
solutions with de Sitter branes. First, we investigate thermodynamic stability in terms of de Sitter 
entropy. We see that the family of solutions is divided into two distinct branches: the high-entropy 
branch and the low-entropy branch. By analogy with ordinary thermodynamics, the high-entropy 
' branch is expected to be stable and the low-entropy branch to be unstable. Next, we investigate 

. dynamical stability by analysing linear perturbations around the solutions. Perturbations are de- 

■ composed into scalar, vector and tensor sectors according to the representation of the 4D de Sitter 

£\j ' symmetry, and each sector is analysed separately. It is found that when the Hubble expansion 

rates on the branes are too large, there appears a tachyonic mode in the scalar sector and the 
background solution becomes dynamically unstable. We show analytically that the onset of the 
thermodynamic instability and that of the dynamical instability exactly coincide. Therefore, the 
, braneworld model provides a new example illustrating close relations between thermodynamic and 

dynamical instability. 

' PACS numbers: 04.20.-q, 04.70.Dy, 04.50.+h, 11.25.Mj 

. I. INTRODUCTION 

o 

The braneworld scenario, in which matter fields are confined on a brane in higher dimensional spacetime, has been 
one of the most exciting subjects in the research into the early universe. As the early universe is supposed to have 
experienced a rather high-energy epoch, studies of the dynamics of the early universe might reveal as yet undiscovered 
imprints of extra dimensions to us. 
Oh, In our previous studies [l|, 0, HJ we considered a 6-dimensional braneworld model with one or two 3-branes. The 
shape and the size of the extra dimensions are stabilized by magnetic flux of a U(l) gauge field and the geometry 
on each brane is either Minkowski or de Sitter. Since all moduli are stabilized, it is expected that the 4-dimensional 
Einstein gravity should be recovered on each brane. In Ref. [l[ we analyzed how the Hubble expansion rate on each 
brane changes when the brane tension changes. The resulting relation between the Hubble expansion rate and the 
brane tension can be considered as an effective Friedmann equation, which in turn defines an effective Newton's 
constant. It was shown that at sufficiently low Hubble expansion rates, the effective Newton's constant resulting from 
this analysis agrees with that inferred by simply integrating extra dimensions out. The stability of this model was 
investigated in Ref. 0, [I[ , where it was shown that the background with Minkowski branes is always stable against 
linear perturbations with axisymmetry in the extra dimensions. 

The purpose of this paper is to study stability of backgrounds with de Sitter branes. This is not just a trivial 
extension of the previous analysis but actually results in a qualitatively different picture. Indeed, the background 
becomes unstable when the Hubble expansion rate becomes too large. (For other models, see for example 0J, MM !•) 
We shall obtain the critical expansion rate at which the background becomes unstable and interpret this behavior in 
terms of de Sitter thermodynamics. 

Spacetimes with horizons appear mysterious in many respects. An observer does not have causal contact with any 
events beyond a horizon and, in this sense, the horizon separates the region behind it from the observer's side. Thus, 
from the observer's viewpoint, information beyond the horizon is completely irrelevant (as far as classical evolution in 
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the observable region is concerned). This, nothing but a restatement of the definition of horizons, suggests a coarse- 
grained description of a physical system in or of a spacetime with horizons. It is therefore not totally surprising that 
a spacetime with horizons exhibits properties analogous to those of thermodynamics. Indeed, in the case of black 
holes such properties are known as black hole thermodynamics . 

Black hole entropy, defined as one quarter of horizon area in Planck units, is central to black hole thermodynamics. 
The second law of black hole thermodynamics states that black hole entropy never decreases classically. Its semi- 
classical extension, known as the generalized second law, states that the sum of black hole entropy and matter entropy 
does not decrease. These second laws are direct analogues of the second law of ordinary thermodynamics and restrict 
directions in which a physical system involving black holes can evolve. In ordinary thermodynamics it is due to the 
second law that entropy can be used to analyze stability of a thermodynamic system. Therefore, the validity of the 
second law of black hole thermodynamics suggests potential use of black hole entrop y to guess stability of a system 
including black holes or objects (branes, strings, rings, and so on) [H HO, El El GSM El- (Also, for a review see 
e -g- [13 an d references therein.) 

The close connection between black hole event horizons and thermodynamics can be extended to event horizons in 
cosmological models with a positive cosmological constant [la |. The de Sitter thermodynamics includes concepts of 
de Sitter temperature and de Sitter entropy. 

In this paper we shall extend the de Sitter entropy to our braneworld set-up with two extra dimensions and use it to 
"predict" the stability condition for the backgrounds with de Sitter branes. Intriguingly, the simple "prediction" by 
de Sitter entropy completely agrees with the result of rather involved analysis of linear perturbations. This is really 
non-trivial and may be considered as an evidence showing the usefulness of de Sitter entropy and, more generally, 
holographic principle [H, H3| • 

The rest of this paper is organized as follows. In Sec. |TT] we briefly review our 6-dimensional braneworld model 
with de Sitter branes. In Sec. IIIII we study thermodynamic stability of the spacetime by using entropy. In Sec. IIV1 in 
order to investigate dynamical stability we consider linear perturbation around the background solution. In Sec. |V| 
we summarize the results and discuss them. 



II. 6D BRANE MODEL WITH WARPED FLUX COMPACTIFICATION 



We consider a 6D Einstein-Maxwell system described by the action 
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where Ag is the 6D bulk cosmological constant (Ag > 0) and Fmn = 9m An — On Am is the field strength of the U(l) 
gauge field Am introduced to stabilize the extra dimensions. (We use units with Gq = 1 where Gq is the 6D Newton's 
constant unless otherwise noted.) The braneworld solution is 

df 2 - 

ds 2 = f 2 ds 2 + — + /(W 2 , A M dx M = A[r)d(j>, (2) 



where 



A 6 ~ 2 M b 2 7 ,_ b 



and ds 2 is the metric of the 4-dimensional spacetime with a constant curvature k, i.e. de Sitter (k > 0), Minkowski 
(k = 0), or anti-de Sitter (AdS) (k < 0) spacetime. Without loss of generality, we can normalize k to ±1 except for 
the Minkowski (k — 0) case. We consider the range of parameters in which /(f) has two positive roots, and locate 
two 3- branes at the positions of the two roots f — r± (0 < r_ < r + ). We call the branes at f = r± the f±-branes, 
respectively. The tension of each brane determines a conical deficit in the extra dimensions. To be precise, the period 
of the angular coordinate (<j) ~ <j) + A(j>) is given by 

2tt - 8.q + = 2tt - 8^x_ 
\f(r + )/2\ |/'(r_)/2| 



For a review see e.g. [i| and references therein. 
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where a± are tensions of the branes at r = r± . 

This brancworld model includes two compact extra dimensions and two 3-branes at r = r± with the tensions a±. 
For k = +1, the geometries on the brane at f — r± are the 4D de Sitter spacetimes with the Hubble expansion rates 
h± = l/r±, respectively. For k = 0, the brane geometry is the Minkowski spacetime. In each case the energy scale 
on each brane is scaled by the warp factor r 2 depending on the position of the brane. The role of the U(l) field is to 
stabilize extra dimensions with magnetic flux. 

Now, for later convenience we re-scale coordinates and parameters of the geometry with the 4D Hubble parameter 
h + on the r + -branc in the following manner: 

r - M. ~ b 

f^-, (f)^h<j), , b -» j-j, x^^hx", (5) 

where we have omitted the subscript "+" of h + and the tilde denotes re-scaled quantities. After this re-scaling, we 
obtain 

dr 2 

ds 2 = r 2 g^dxW + J^+ f{rW, A M dx M = A{r)&4>, (6) 

where 

Here g^ v is the metric of the 4D de Sitter spacetime with the Hubble parameter h. The metric function now satisfies 
f(a) = /(l) = 0, where a = r-/r+ is the ratio of the warp factors at the branes and, at the same time, describes how 
the geometry of the 2D extra dimensions is warped. Note that in the special case a = 1 the geometry is not warped 
but locally a round sphere. 

Advantages of this re-scaled parameterization are as follows. First, the new parameters and coordinates remain 
finite and regular in the limit where the geometries on the branes become Minkowski. Second, since the warp factor 
evaluated at the r + -brane is unity (r = 1), many quantities rooted in these coordinates have physical meaning as 
observables measured by observers on the r + -brane. As we will see later, a Kaluza-Klein (KK) mass measured on the 
r + -brane remains finite in the a — > limit, while that on the r_-branc diverges. Hence, it is useful (in particular for 
the purpose of numerical calculations) to use quantities measured on the r + -branc. Of course, if we know the value 
of a quantity measured on the r + brane then we can easily obtain the corresponding value measured on the r_ -brane 
by multiplying the warp factor a, e.g., h- = h + /a, m 2 _ = m+/a 2 7 and so on, where m± is the KK mass observed on 
the r±-brane. 

We introduce another convenient parameter /3. From f(a) = /(l) = we obtain the relation 

2 b 2 _ A 6 (q 4 + a 3 + a 2 + a + 1) 

12a 3 ~ 10(a 2 + a + l) ' ( > 

Since b 2 > 0, the Hubble parameter h has a maximum value h max (a) for a given a. Hence we define /3 by 

a- h " _ 10(1 + a + a 2 ) h2 (Q) 

/imax(a) 2 A 6 (l + a + a 2 +a 3 + a 4 ) ' 

which is the squared Hubble parameter h 2 on the r+-brane normalized by h max (a). Alternatively, (3 could be defined 
as h 2 _ normalized by its maximum value, but the two definitions are actually equivalent. Thus, (3 describes how much 
the 4D part of the metric is curved. When (3 = 0, the 4D spacetime becomes flat 2 . In contrast, in the maximal case 
((3 = 1), the flux disappears: 6 = 0. 

For a fixed bulk cosmological constant A 6 and a fixed discrete parameter k (= 0, 1), the solution is locally parame- 
terized by the two parameters b and M.. However it is more convenient to use the pair of dimension- less parameters 
(a,/3), both of which run over the finite interval [0, 1]. This parameterization includes both k = 1 (0 < (3 < 1) and 
k = (fl = 0) cases, and treats the de Sitter and the Minkowski branes in the common way. 

In the following sections we use the parameters (a, (3) and the coordinates (r, 0, x M ) to describe the background 
geometry. 



2 In fact, if (3 < then the solution represents a braneworld with AdS branes. In this paper, however, we shall not consider this case, 
focusing only on the de Sitter (0 < fi < 1) and Minkowski (/3 = 0) branes. 
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III. THERMODYNAMIC STABILITY 



As stated in the introduction, one of the main subjects in this paper is the close connection between thermodynamic 
and dynamical properties of spacetimes. In particular, we shall see exact agreement between thermodynamic and 
dynamical stabilities for the 6-dimensional brane world solutions described in the previous section. In this section 
we discuss thermodynamic stability by using the de Sitter entropy. We also point out the close relation between the 
specific heat and the effective Friedmann equation. 



A. Entropy argument 

In the 4-dimensional Einstein gravity, a de Sitter space has entropy given by one quarter of the area of the hori- 
zon [l8| . This is an analogue of the well-known black hole entropy and has been the basis of holographic arguments 
for de Sitter spaces [H, 0] ■ 

In the 6D braneworld background solution, each point in the 2D extra dimensions corresponds to a 4D de Sitter 
space. Therefore, we shall define the de Sitter entropy in our context by one quarter of the area of the cosmological 
horizon integrated over the extra dimensions, i.e., one quarter of the volume of the 4-dimensional surface foliated by 
de Sitter horizons: 

where A4 is the 4-dimcnsional volume given by 

4-K f 1 /" A * 47T 

M = ^ J r 2 dr J ^=—(1- «3 )A ^ (11) 

and Gq is the 6D Newton's constant. Note that, for the convenience of the discussion, we have temporarily restored 
Gq in this paragraph. Actually, it is easy to show that this definition agrees with the definition of the 4D de Sitter 
entropy that observers on each brane would adopt. The area of the de Sitter horizon and the Newton's constant on 
the r± -brane are, respectively, given by 

A+ = a- 2 A- = % (12) 



h 2 



and 



1 a 2 1 



2 1 r i 



Gn+ Gn- Gq J a 
Thus, the 4D de Sitter entropies on the r±-branes are 



— / r 2 dr / d<f>. (13) 



S± = (14) 
4Gn± 

respectively. These do indeed agree with the above definition based on the horizon area integrated over the extra 
dimensions: 

S+ = S-= S. (15) 

In the 4D picture, the integration over extra dimensions is taken care of by the above formula of the 4D Newton's 
constant (|13p . Since the three different definitions agree, we have a unique definition of de Sitter entropy in our 
braneworld set-up. 

To investigate the thermodynamic stability of the system, it is useful to consider the de Sitter entropy as a function 
of conserved quantities. The family of solutions described in the previous section has three conserved quantities: the 
total magnetic flux and the tensions of two branes. The total magnetic flux $ is given by 

*= f F riP dr I d0=^(a 3 -l)A0. (16) 
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By Eq. ([4}, tensions a± of the r±-branes are expressed as 



respectively, where n± are defined by 



2tt-k ± 



K+ = ~f , (l)A4>, «_ = i/'(a)Af (18) 



Since there is one-to-one correspondence between <r± and k±, we can consider the de Sitter entropy S as a function 
of either ($, cr + , <r_) or ($, k + , k_). 

Variables S, $ and /cj- defined above are proportional to A<f>. Because of this trivial dependence, we can eliminate 
A(f) from all thermodynamic considerations by properly normalizing the variables. In particular, in Appendix [X] it is 
shown that the entropy normalized by K-, S = S/k-, satisfies the differential relation 

M=(^)^+(£)d$, (19) 



V3/i 4 / ' \3h 4 , 

where r\ = and $ = is the normalized total flux. This means that the normalized entropy S is a function 

of the two variables r\ and $: S(r), $). Hereafter, r\ and $ are to be considered as two independent thermodynamic 
variables. This accords with the observation made in the previous section: the family of solutions are parameterized 
by two parameters (a, 0) after proper scaling of variables. The angular period A(f> has been eliminated so that S, r\ 
and $ do not depend on it. 

Note that the normalization factor 1/k— is just a fixed constant from the viewpoint of observers on the r + -brane 
since these observers are to see the response of the 4D geometry induced on the r+-brane to the change of the tension 
<T+ or equivalently, k+. In particular, with k_ fixed, we have 

dn = - — da + . (20) 



Thus, the relation (|19p is relevant for observers on the r+-brane (but not for those on the r_-brane). Alternatively, 
we could write down a differential relation relevant for observers on the r_-brane by normalizing the variables by 
the factor 1/k+, which is a constant from the viewpoint of observers on the r_-brane. In this sense, we have two 
conceptually different pictures: one from r+-brane perspectives and the other from r^-brane perspectives. The result 
of the thermodynamic stability analysis does not depend on the picture: the result from one picture completely 
agrees with that from the other picture. In the following, we shall see the thermodynamic stability from the r + -brane 
perspectives. Thus, unless there is a possibility of confusion, we omit the subscripts "±" . 

Fig. [1] shows the graph of the normalized de Sitter entropy S as a function of the conserved quantities r\ and 3>. 
As we can easily see, S(r], $) is not single-valued but actually double-valued in some part of its domain. It is also 
easy to see that derivatives of S(rj, $) are singular on a part of the boundary of the domain. We call this part of the 
boundary of the domain the critical curve. Different points on the surface S = S(rj, $) represent different background 
solutions, i.e. different spacetimes. (There is no background solution outside the domain of S(rj, $).) Therefore, the 
critical curve divides the family of solutions into two distinct branches: one with larger S and the other with smaller 
S. We call them the high-entropy branch and the low-entropy branch, respectively. By definition, these two branches 
merge on the critical curve. Note that in some region of the low-entropy branch, S(r), $) is single- valued. 

If the de Sitter entropy has any physical meaning like entropy in ordinary thermodynamics then the high-entropy 
branch must be thermodynamically preferred in some sense and must be more stable than the low-entropy branch. 
The spacetime is marginally stable on the critical curve. These observations are based only on our intuition that the 
de Sitter entropy should play a role similar to that of entropy in ordinary thermodynamics. Nonetheless, surprisingly, 
we shall see in the next section that these observations correctly "predict" the result of rather detailed analysis of 
dynamical stability against linear perturbations. 

We now derive an explicit equation describing the critical curve. For this purpose, we consider the coordinate 
transformation from (a, (3) defined in the previous section to (77, $) defined above. As we shall see below, the critical 
curve is the set of points where this coordinate transformation becomes singular. As noticed in Sec. HH a (s (0, l])is 
the ratio of the warp factors of two branes and describes the shape of the two-dimensional extra space. On the other 
hand f3 (€ [0, 1]) represents the curvature of 4D de Sitter space normalized by the maximum value ft^ ax (a) for a given 
a. Since the background solution is uniquely specified by (a, (3) up to the scaling explained in the previous section, 
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FIG. 1: The normalized entropy S as a function of r\ and Different points on the surface S — S(rj, $) represent different 
background solutions, and there is no background solution outside the domain of S(rj,<&). It is easily seen that S(rj, $) is not 
single- valued but actually double- valued in some part of its domain. Thus, the family of background solutions has two branches, 
divided by a critical curve on which some derivatives of S(ri, $) are singular. The high-entropy branch is thermodynamically 
preferred to the low-entropy branch. 



any (properly normalized) quantities can be rewritten as functions of (a,/3). In particular, we have two single- valued 
functions rj(a, f3) and <E>(a,/3) as a map (a, (3) — > (77, <E>) . The entropy S is also a single- valued function of (a,/3) and, 
thus, the critical curve is determined by vanishing Jacobian determinant of this map: 



9(77,$) 



= 0. 



This condition can be solved with respect to (3 to give the following expression for the critical value of (3: 

P CTi (a) = Wl , 1 1 „ 2 , 2/1 1 1 , „ 3 [11 + 33a + 66a 2 + 85a 3 + 90a 4 + 85a 5 + 66a 6 + 33a 7 + 11a 8 



(21) 



12(1 + a + a 2 ) 2 (l + a + a 2 + a 3 4 - , 
- (1 + a)(l + 2a + 4a 2 + 2a 3 + a 4 ) \/\ + 34a 3 + a 6 ] . 



(22) 



This is the analytic expression for the critical curve. As we see later, /3 cri (l) = 2/3 (or, h 2 ri = A 6 /9 and b 2 ri = 6/i 2 ri ) 
is nothing but the critical value for stabilizing the extra dimensions like in Freund-Rubin compactification Q ■ 



B. Effective Newton's constant 



We have seen that the (normalized) de Sitter entropy is not a single-valued function but actually a double-valued 
function of the conserved quantities (ry, $) in some region of the domain. This observation has led us to a rather 
natural stability criterion that the high-entropy branch should be stable and that the low-entropy branch should be 
unstable. This is a global statement: we would not be able to say one of the two branches should be stable or unstable 
without knowing the existence of the other branch. 

In this subsection, we look at the stability from a different viewpoint, considering the effective Friedmann equation 
investigated in Ref. [l[. In particular, we shall see that the thermodynamic stability is equivalent to the positivity of 
an effective Newton's constant. Note that, while the thermodynamic stability is a global statement, the positivity of 
the effective Newton's constant is a local statement. 

In our previous work [l[ we showed that the 4D Friedmann equation is recovered on the brane at sufficiently low 
energy as the response of the 4D Hubble expansion rate to the change of the brane tension. To be more precise, for 
small Hubble expansion rates h, h 2 is expanded as 



[a-a ) + O{{a-(j Q f), 

a=ao 



(23) 



FIG. 2: The Hubble expansion rate squared, h 2 , as a function of rj and $. For small h 2 , (dh 2 / '9rj)g is negative. However, as 
h 2 increases, the slope turns over and (dh 2 / drf)^, becomes positive. (This figure is identical to Fig. 3 of Ref. Q.) 



where a is the brane tension and the 4D Newton's constant Gn is obtained by simply integrating extra dimensions 
out as (fl3|) . The constant oo is the value of a for the Minkowski brane (h = 0) and depends on $ and er_ , which are 
actually constants from the viewpoint of observers on the r + -brane. Note that the 4D Newton's constant Gn defined 
as (1131) is always positive. 

The terms of order 0((a — oo) 2 ) an d higher are small compared with the linear term if the Hubble expansion rate 
is sufficiently lower than the energy scale set by the bulk cosmological constant. However, for larger h 2 , the higher 
order corrections become relevant and h 2 is no longer approximately linear. Moreover, as we can see from Fig. [2| h 2 
is not necessarily an increasing function of the brane tension. Indeed, if h 2 is larger than a critical value then it is a 
decreasing function of the brane tension. To make this peculiar behavior more quantitative, let us define the effective 
Newton's constant G e ff on the r + -brane by 



and characterize the critical value by 1/G e ff = 0. For h smaller (or larger) than the critical value, G c g is positive (or 
negative, respectively). 

Note that, since this effective Newton's constant runs as the Hubble expansion rate h changes, the exact agreement 
with Gn is met only in the h 2 — > +0 limit. The agreement in this limit itself is a rather non-trivial consequence 
of the dynamics of bulk spacetime as explicitly shown in 1], but physical interpretation is simple. While Gn was 
determined via the coefficient of the contribution of the graviton zero mode to the higher dimensional action, G c g here 
incorporates contributions of all Kaluza-Klein modes as well as the zero mode. Since Kaluza-Klein contributions are 
irrelevant at low energies, these two definitions should agree. On the other hand, at high energies the Kaluza-Klein 
modes become relevant and make G g deviate from Gn- 

In order to examine the curve 1/G e ff = 0, the following observation is useful: 




(24) 




6h 

«_ d(a,(3) 



I 



d(a,P) 



(25) 



Since h and $ are single- valued smooth functions of (a, /?), the Jacobian \d(h, $)/9(a, (3) \ does not diverge. Hence, G c g 
diverges if and only if the denominator \d(r], <£>) / d(a, f3) \ vanishes. This condition is nothing but Eq. (|2Tj) characterizing 
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the critical curve, which divides the parameter space of background solutions into the high-entropy branch and the 
low-entropy branch. Therefore, 

Critical curve <S=> 1/G c s — 0, 
High-entropy branch l/G e g_ > 0, 

Low-entropy branch 1/G c s<0. (26) 

Additionally, it is worth mentioning that the effective Newtion's constant has a connection with thermodynamic 
quantities on the brane. In the 4D case we can define a temperature of the de Sitter space by the Hubble parameter 
H as T = H/2n. Then, in a straightforward manner, the specific heat is given by 

dS\ _ 2tt 



T ll^) =~T^m- (27) 



Now, concerning the 4D de Sitter brane, we consider a specific heat regarding T — h/2ir as its temperature. From 
Eq. 11^1) we obtain 

rr,f9S\ , / dS\ ms_ /0n\ 2tt 

T L = K -M^ =w«r= -77-15- (28) 



dTj $a _ \ dh )$ 3/l3 \dhJs G cS h 2 



This expression is the same as that of the 4D de Sitter space except for its gravitational "constant." For the 4D de 
Sitter brane the sign of the specific heat will change on the critical curve also. Hence, from the viewpoint of the 4D 
effective theory, we conclude that the effective Newton's constant on the brane determines thermodynamic properties 
of the braneworld. 



IV. DYNAMICAL STABILITY 



In the previous section we have discussed the thermodynamic stability of the 6D braneworld model as a property 
of the background solution. Now let us investigate the dynamical stability of this spacetime. For this purpose, we 
consider linear perturbations around the background solution. Perturbations are decomposed into scalar-, vector- and 
tensor-sectors according to the representation of the 4D de Sitter symmetry. If the mass squared of the Kaluza-Klein 
modes m 2 is negative, then the corresponding mode is tachyonic and unstable. 

Since the background spacetime forms a two-parameter family of solutions, the system of the perturbation equations 
depends on two parameters (a,/3). We investigate the lowest mass squared m 2 (a,/3) in each sector as a function of 
the two parameters in the square domain 0<a<l,0</3<l. Our strategy for attacking the problem is as 
follows. First, we analytically solve the perturbation equations in the a = 1 case for arbitrary [3 € [0, 1], and obtain 
m 2 (a = 1,13) on the side a = 1 of the square domain. Next, by successively applying the relaxation method, we 
numerically calculate m 2 {a,l3) as a changes from 1 to for every f3. In each step of the second procedure, the result 
of the previous step with a slightly larger value of a is used as an initial guess for the relaxation method. 



A. Scalar perturbation 

As shown in Appendix [C] the scalar perturbation with an appropriate gauge choice is given by 

dr 2 

g M Ndx M dx N = r 2 (l + ^ 2 Y) 9liV dx^dx v + [1 + ($! + $ 2 )Y]— + [1 - ($ x + 3$ 2 )Y}fdtf, 

J (29) 

A M dx M = (A + a F)d</>, 

where Y is the scalar harmonics on the 4D de Sitter space satisfying 

V 2 y - m 2 Y = 0. (30) 

Here, V M is the covariant derivative associated with the 4-dimensional de Sitter spacetime with the Hubble h. The 
Einstein equation and the Maxwell equation are reduced to the following two perturbation equations for <I>i and $2: 



„lf ^ 4A 6 ^ ^ , m 2 + 18/i 2 
+2 ( ±j + - *i - + $2) + — f $1 = 0, 

4 . m 

^ 2 + 27V 



( 31 ) 



$2 + r $ 2 + 7^77(^1 + 2*2) = 0, 
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where the prime denotes the derivative with respect to r, and m 2 is the eigenvalue of the harmonics. 

The boundary conditions at r = a, 1 are obtained by setting the coefficients of 1// in (j3"Tj) to zero at the brane 
positions, as in (2j. Hence we have 



2/'$; - 4A 6 ($i + $ 2 ) + 



$1 



($i +2<J> 2 ) 



0. 



(32) 



r— a,l 



An alternative and more rigorous derivation of the boundary conditions is to use the formalism developed in [l| . The 
result is 



/*llr=a,l = [2/ / *a + (/#l) / ] P=Q ,i = 0. 



(33) 



These two sets of boundary conditions lead to the same Taylor expansion of $i and $ 2 in the neighborhood of the 
boundaries and, thus, are equivalent. 

In the original coordinate system (r, <fi) the geometry in the a = 1 limit appears to be singular because two positive 
roots of f(r) corresponding to the brane positions coincide. Actually, the proper distance between the r±-branes 
remains finite and it turns out that the apparent singularity due to coordinate artifacts is removed by appropriate 
coordinate transformations. Indeed, in the coordinate system (w,ip) defined by 



2r - (1 + a) 

w = , ip={l-a) 

1 — a 



(34) 



the geometry is obviously regular in the a = 1 limit. Note that the coordinate w ranges over the interval [—1, 1]. 
Since the branes are always located at w = ±1 for any a (even for a = 1), the new coordinate w is more useful for 
numerical calculation than the original one r. Therefore we rewrite the perturbation equation in terms of w as 



5(1 - a) 



d w _f _ 
f [(l-a)w + (l + a)} 
02^ 4(1 -a) 



[(1 - a)w + (1 + a)} 



d w <S>i - j($i 
d w <$>2 + 
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2/i 2 + I8h 2 /A 6 



[(1 - a)w + (1 + a)] 2 / 
,2 

? ($! + 2$ 2 ) = 0, 



(35) 



[(l-a) W + (l + a)] 2 / v 



where 



A 6 (l-a) 2 
1 

"40(1 - a) 2 



-[{1 - a)w + {1 + a)} 2 + 



2 4(l + a + a 2 +a 3 + a 4 )/3 



1 + a + a 2 



(36) 



32(1 + a)[l + a 2 (l + a 2 + a 4 )(l - 0) - (1 + a 3 
[(1 - a)w + (1 + a)] 3 (l + a + a 2 ) 



256a 3 (l + a + a 2 + a 3 + a 4 )(l - /?) 
[(1 - a)w + (1 + a)] 6 (l + a + a 2 ) 



and 



m 

2K' 



The boundary condition is written as 

2d w fd w $ 1 - ($! + $ 2 ) 



r$i 



2/i 2 + 18fe 2 /A 6 
[(1 - a)w + (1 + a)} 

($ x +2$ 2 )| w=±1 



0. 



w=±l 



= 0. 



(37) 



(38) 



1. Analytic solution for a = 1 



Now we show the analytic solution for a 
a round 2-sphere as in Ref. 



1. In this case the geometry of the extra dimensions becomes locally 
2l| and the warp factor is a constant. Hence, the background spacetime for a = 1 
corresponds to the football-shaped extra dimensions [13, [HI . 
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In the a = 1 limit, the perturbation equations (|35[) are reduced to 



(1 - - W a - (2 - ^L_) * - = 0, 

(1 - w 2 )d 2 w $ 2 + + 2$ 2 ) = 0. 



(39) 



4-3/3 

and the boundary conditions are given by 

-2w(4-3/3)$' 1 -4$ 2 + (-4 + 3/3 + 2 At 2 )$ 1 | iu=±1 = 0, $i + 2$ 2 | M=±1 = 0. (40) 
The solution satisfying these boundary conditions is 

^ =-2P n (w) + ^^-[n(n+l)P n (w)-2wP^(w)], $ 2 ±} = P n (w), (41) 

up to an overall constant, where n = 0, 1, 2, • • • , and P n is the Legendre function of the first kind. We obtain the mass 
spectrum of the scalar perturbation as 



2 _ 2 _ 2 - 3/3 + n(l + n)(4 - 3/3) ± y/(2 - 3/3) 2 + 12n(l + n)(l - gjg - 3/3) 
M - A*(±) = 1 ■ (4^J 

The mass of the lowest mode is given by 

M 2 +) (n = 0) = (43) 
Thus, for /3 > 2/3, the lowest /Lt 2 becomes negative and the background spacetime is destabilized!!, [f| 0]. 



2. T/ie lowest mass of the KK mode 



Now we numerically solve the perturbation equations and obtain the lowest KK mass of the scalar mode for general 
a. Using the analytic solution for a = 1 as a first initial guess, we solve the problem for a slightly smaller value 
of a by the relaxation method. Then, as we change a to smaller values towards step by step, we in turn use the 
numerical solution of the previous step as an initial guess for the relaxation method and obtain the mass squared for 
each a and /3. 

Fig. [3] shows m 2 for the lowest mode of the scalar perturbation. For (3 = 0, which means the 4D spacetime is the 
Minkowski, we have reproduced the result of our previous study Q that there is no unstable (i.e. neither massless nor 
tachyonic) mode. We find that as /? becomes larger, namely as the Hubble parameter on the brane becomes larger, 
in 2 decreases and eventually becomes negative. Therefore we conclude that when the Hubble on the brane is too large 
the extra dimensions are destabilized and such configurations are unstable for general a. In addition, in the case of 
vanishing flux ((3 = 1) the configurations are always unstable. 



3. The threshold massless mode 



Supported by the result shown in Fig. [3l the assumption that m 2 is a continuous function of a and (3 is made. 
Hence, it must vanish on the boundary between the stable and unstable regions. In other words, there should appear 
a massless mode as one moves from the stable/unstable region to the unstable/stable region. 

We can actually obtain this threshold massless mode analytically. The general solution of (|3"Tj) with m 2 = is given 
by 

5E + 



I (2L -l Si 

f Vr 3 r 6 



c 3 r 



£4 
r 



— , $2 



(44) 



where Ci(i = 1,2,3,4) are arbitrary constants. By imposing the boundary conditions (|33p . we will be able to obtain 
the boundary between the stable and the unstable region on the (a, f3) plane, where the massless modes appear. The 
boundary conditions lead to 



a 2 


a- 1 


a 7 




a 4 




1 


1 




1 


6A 6 a 6 


12A 6 a 3 


-206 2 a 3 + 4A 6 a n 


- 120A^a 6 


-56 2 + 4A 6 a 8 


6A 6 


12A 6 


-206 2 + 4A 6 - 


120M 


-5b 2 + 4A 6 



30 Ma 3 
30M 










) 








c 3 






w 





(45) 
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FIG. 3: m 2 for the lowest mode of the scalar perturbation. When (3 = the 4D spacetime is the Minkowski. As /3 becomes 
larger (i.e. the 4D Hubble becomes larger), m 2 decreases and eventually becomes negative beyond some value of f3. The line 
represents the boundary between the stable and the unstable region, where the massless modes appear. 



where M. and b are functions of (a, (3) determined by f(a) = f(l) = 0. In order for this equation to allow a non-trivial 
solution for (01,02,03,04), the determinant of the matrix must vanish. From this condition, we obtain the relation 
between (a,/3) as 



B = (3 cli (a) = —— n , , 1 - - — 111 + 33a + 66a 2 + 85a 3 + 90c/ + 85a 5 + 66a 6 + 33a 7 + 11a 

' 12(1 + a + a 2 ) 2 (l + a + a 2 + a 3 + a 4 ) L 



s 



- (1 + a)(l + 2a + 4a 2 + 2a 3 + a 4 ) \J\ + 34a 3 + a 6 ] . 

(46) 

This relation is exactly identical with (|22j) , which defines the critical curve on which the high-entropy branch and the 
low-entropy branch merge. The threshold massless mode divides the whole square (a, /3)-domain into positive and 
negative m 2 regions, while the critical curve in Sec. Mil divides the whole square (a, /3)-domain into high- and low- 
entropy regions. Therefore, the high-entropy and the low-entropy branches in Fig. [1] exactly agree with the positive 
and negative m 2 regions in Fig. [3l respectively. This fact implies that in this spacetime thermodynamically preferred 
configurations are dynamically stable, while thermodynamically non-preferred ones have tachyonic modes and become 
dynamically unstable. The onset of dynamical instability coincides with the onset of thermodynamic instability. This 
is similar to the Gubser-Mitra conjecture for extended black objects fl3l. IT3|. 



B. Vector perturbation 

As shown in Appendix [Cl the vector perturbation with an appropriate gauge choice is given by 

d^ 2 

g M Ndx M dx N = r 2 g^dx^dx» + 2/ l(T)0 V (T)M d^d0 + — + /d</» 2 , 

A M dx M = a (T) V (T)M da; AI + Ad<f>, 

where Vrrw is the vector harmonics on the 4D de Sitter space, which satisfies the transverse condition, VVrrw = 0. 
The Einstein equation and the Maxwell equation are reduced to two perturbation equations with respect to /i(T)^ an d 
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0(T), 



26 , 



6h \ n 



r 2 / 



(r 2 /a' (T) )' + 6 



l (T)<t> \ 
r 2 J 



+ m 0( T ) = 0, 



where m 2 is the eigenvalue of vector harmonics on the de Sitter space with the Hubble constant h: 



VV (T) , - (m 2 + 3h 2 )V {T) ^ = 0, VV (t)a 



0. 



(48) 



(49) 



It is again convenient to use the coordinates (w,ip) defined in Eq. (|34[) which do not have any singular behaver in 
the limit a — > 1 . The perturbation equations and the boundary conditions requiring regularity at the brane positions 
are written in the coordinates (w,ip) as 



1 



r 6 d u 



1 



M 2 + 3h 2 /A 6 



2r 2 f 



4-9 



d w {r 2 fd w * l ) + h -d w 



(50) 



and 



( M 2 + 3/i 2 /A 6 )* 2 | 



r 2 d w fd w *i + -3 U 



= 0, 



tu=±l 



where 



*1 



6a 



(T), 



n (T)<P 

1-a' 



(51) 



(52) 



Note that the factor (1 — a) -1 in the definition of iff 2 reflects the fact that the relation between the original angular 
coordinate </> and the new coordinate tp includes the factor (1 — a): h^^dcj) = ^dy. On the other hand, the factor 
6 in the definition of \&i is not essential but is just to absorb the factor 6 in the perturbation equations. 



1. Analytic solution for a = 1 
For a = 1 the perturbation equations become 

04(1 ~ ^ 2 )d^i] + 8 , (1 ~ w d w *2 + t^L^ i 0- 



4-3/5 Jw 2 • 4-3/? 



(l~w 2 ){dl* 2 -d w * 1 } + 



4/z 2 + 2/3 



#2 = 0, 



4-3/3 

where /i 2 = m 2 /2A 6 . The solution satisfying the boundary conditions is given by 

^ ±} H = (/3 ± V/? 2 + 8n(l + n)(l - /3)(4 - 3/3)) P n (tu), 

^(±) M = (4-3/3)(l-^ 2 )P T ;H 
2 /3 ± ^/3 2 + 8n(l + n)(l - /3)(4 - 3/3) ' 

up to an overall constant, and the mass spectrum of the vector perturbation is 



2 _ 2 _ n(l + n){A - 3/3) - /3 ± V/3 2 + 8n(l + n)(l - /3)(4 - 3/3) 
M - M( ±) = t 



(53) 



(54) 



(55) 



where n = 0, 1,2, • • • for fi 2 + ^ and n = 1,2, • •• for A* 2 _)- There are two lowest mass modes for /i 2 + ^(n = 0) and 
li 2 _^(n — 1), and both are massless. In other words the vector perturbations have two zero modes: 

4-3/3, 



= 6d( T ) = Cl + 8c 2 (l - /3) W, * 2 = fr(T)v 



C2" 



-(1-^), 



(56) 



where C\ and c 2 are arbitrary constants. As shown below for general a, these two zero modes represent physical 
degrees of freedom of the Maxwell field and the gravi-photon. 
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FIG. 4: m 2 for the first two KK modes of the vector perturbation. When a = 1, the values correspond to /i 2 + j for n = 1 (left) 
and /i(N for n = 2 (right). These are non-negative for the entire region of (a, (3). 



2. Zero modes for general a 

For general a we find two zero modes corresponding to the above ones as 

O(T) = Ci + C 2 A /l(T)0 = c 2 /, (57) 

where ci and c 2 are arbitrary constants. It is obvious that the degree of freedom represented by c\ originates from 
the original 6D Maxwell field and, thus, is associated with the 4D part of the U(l) gauge symmetry. On the other- 
hand, as seen below, the other degree of freedom represented by c 2 is associated with coordinate transformation of 
the angular coordinate <f> and, thus, can be regarded as a gravi-photon. 

We now see this explicitly. Let us consider an infinitesimal U(l) gauge transformation represented by a parameter 
x(x) and an infinitesimal coordinate transformation cj> — > <f> + £(x), where x( x ) an d £( x ) depend only on the 4D 
coordinates. Under these, the /x-components of the perturbation of the U(l) gauge potential and the (/z</>)-components 
of the metric perturbation transform as 

SAp -> SAp + d^ x + Ad^, -> + fd M £, (58) 
and all other components are unchanged. Thus, if we define two 4-vectors and as 

6A» = aM+AaW, V = /< 2) , (59) 
then the above transformation law is rewritten as two separate U(l) gauge transformations: 

a^^a^+d^x, a^aP)+8,f. (60) 

Clearly, c\ and c 2 above are coefficients of the transverse components of afi and , respectively. 
Thus, any instability will not occur for these two modes. 

3. Stability and the 1st KK mode for general a 

Having seen that the two zero modes of the vector perturbation remain massless for any a, we will now examine 
the 1st KK modes for general a numerically. Fig. 0] shows m 2 for the first two KK modes of the vector perturbation. 
We find that m 2 is non-negative for the entire region of (a, ft). Therefore, there are no unstable modes in the vector 
sector. 
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C. Tensor perturbation 

The tensor perturbation is given by 



dr 2 

g MN dx M dx N = r 2 (g^ + h (T) T (T) ^dx^dx" + ~+ (gl) 
A M dx M = Ad(p, 

where T(xw is the tensor harmonics on the 4D dc Sitter space, which satisfies the transverse and traceless conditions, 
VT (T)jlu , = T (T) % = 0, and 

V 2 T (t)ai „ - (m 2 + 2/i 2 )T fI> „ = 0. (62) 
The perturbed Einstein equation becomes 

1 



^(r 4 fh{ T) Y + m 2 h {T) =0, (63) 



(64) 



and there is no relevant equation coming from the Maxwell equation. With the coordinate w, this is written as 

- 1 d w {[(1 - a)w + (1 + a)] 4 fd w h (T) ] + 2 M 2 /x (T) = 0. 

[(1 — a)w + (1 + a)J^ 

_Z. Analytic solution for a — 1 
For a = 1 the perturbation equation becomes 

4« 2 

9 w [(l-u; 2 )5 w V)] + 4ZTgV) = °- ( 65 ) 

The regular solution is given by 

h (T) =P n (w), (66) 
where P n is the Legendre function of the first kind. The KK mass squared n 2 is 

9 4-3/3 , 

^ 2 = — i±n(n + l), (67) 
where rt = 0, 1, 2, • • • . It is obvious that the lowest-mass mode is n — and it is massless, representing the 4D graviton. 

2. General a 

For the tensor perturbation it is clear that the zero mode is hm = const., i.e., a homogeneous mode. Fig. [5] shows 
m 2 for the first KK mode of the tensor perturbation. We find that m 2 for the KK modes is non-negative 3 for the 
entire region of (a, 0). Therefore, there is no unstable mode in the tensor sector. (See e.g. (27|.) 

V. SUMMARY AND DISCUSSION 

In this paper we have studied stability of the 6D braneworld solution with 4D de Sitter branes from two different 
perspectives. 



It is known that for tensor modes on the de Sitter background the region m 2 < 2h 2 is forbidden by unitarity 124. |25| . Also, the masses 
in the region 2h 2 < m 2 < ^h 2 are called complementary series |2(I . Now, the masses for the KK modes satisfy m 2 > jh 2 . 
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FIG. 5: m 2 for the first KK mode of the tensor perturbation. The lower plane describes the particular mass, m 2 = 9/i 2 /4, for 
the tensor mode on the de Sitter space. 



One is thermodynamic stability of the braneworld solution. We have defined the de Sitter entropy in the six- 
dimensional braneworld by one quarter of the area of the cosmological horizon integrated over the extra dimensions. 
We have seen that this definition agrees with the 4-dimensional de Sitter entropy defined on each brane. As shown in 
Fig. [1] the de Sitter entropy as a function of the conserved quantities rj and $ is not single-valued but double- valued. 
Therefore, the de Sitter entropy divides the family of solutions into two branches, one with higher entropy (the high- 
entropy branch) and the other with lower entropy (the low-entropy branch), and defines the boundary between them 
(the critical curve). 

The other is dynamical stability of the background solutions against linear perturbations. Perturbations are de- 
composed into the scalar-, vector- and tensor-sectors, according to the 4-dimensional de Sitter symmetry, and we have 
calculated the lowest mass squared in each sector. We have found that when the Hubble expansion rate on the brane 
is larger than a critical value, there is a tachyonic mode in the scalar sector and thus the background is unstable. 
At the critical value, there appears a threshold massless mode in the scalar sector. On the other hand, there is no 
unstable mode in the vector- and tensor-sectors. 

We have found that the critical value at which the threshold massless mode appears is exactly on the critical curve 
dividing the family of solutions into the high-entropy branch and the low-entropy branch. Therefore, we have shown 
that the low-entropy branch is dynamically unstable while the high-entropy branch is dynamically stable. Moreover, 
we have also seen that the 4-dimensional effective Newton's constant is positive in the high-entropy branch and 
negative in the low-entropy branch. In summary, we have shown the equivalence of the following three conditions: 

• Thermodynamic stability. 

• Dynamical stability. 

• Positivity of the effective Newton's constant. 

The close connection between thermodynamic stability and dynamical stability has already been pointed out in 
the literature [l(| EH EH- The result of the present paper, thus, adds yet another example to the list of systems 
exhibiting the close connection between thermodynamic and dynamical properties. In the previous examples in the 
literature, however, one has to rely on numerical analysis to show the equivalence and the boundary between stable 
and unstable regions is not obtained rigorously. On the other hand, our system was simple enough to obtain the 
boundary analytically. In this sense, we may say that the analysis in the present paper is the first example in which 
the equivalence between thermodynamic stability and dynamical stability was proved. Moreover, we have shown that 
for our gravitational system the thermodynamic and dynamical properties of the spacetime is closely related to the 
sign of the effective Newton's constant. 

Non-linear evolution of the unstable solutions is one of the important subjects for the near future. While we 
have shown that the low-entropy branch is unstable, our analysis is limited by the linearized approximation and we 
do not know to what final state the spacetime will evolve from the low-entropy branch. A naive guess is that the 
system should evolve to the corresponding solution in the high-entropy branch with the same values of the conserved 
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quantities, but there remains the possibility that the system might evolve towards a big crunch singularity or evolve to 
another unknown configuration [8|. Indeed, as seen from Fig. [TJ for a solution in the low-entropy branch with a very 
large Hubble parameter close to the maximum value, there is no corresponding solution in the high-entropy branch. 
However, one can show that there is a solution with AdS branes for the same conserved quantities rj and This 
indicates that the FRW cosmology on the brane in this regime should evolve towards a big crunch singularity unless 
there is another unknown stable solution without 4D maximal symmetry. It is worthwhile seeking such solutions 
without 4D maximal symmetry and also analyzing properties of solutions with AdS branes 28]. The whole phase 
structure including the high- and low-entropy branches, the AdS branch and possibly other new branches is unexplored. 

Inclusion of general matter on the brane is also an interesting subject as a future work. In the present set-up, 
the 4D spacetime on the branes has the maximal symmetry and the stress energy tensor on the brane is restricted 
to the form of vacuum energy, or brane tension. It is certainly interesting to investigate dynamical stability and 
thermodynamic properties of the braneworld with more general Friedmann universe and matter contents on the 
brane. For t his purp ose, we need some extensions of the model to include arbitrary matter on the codimension-2 
brane [29l. l30l. l3lTl32lj. It is also interesting to consider supergravity extensions. 

We expect that the close connection between thermodynamic stability and dynamical stability could be extended 
to models with branes of higher codimensions, i.e. models with more extra dimensions. We hope to investigate this 
subject in the future. Here, as the first step, let us briefly discuss extensions to models with more extra dimensions 
but without branes. We consider general Freund-Rubin flux compactifications with dS p x S q . The (p + ^-dimensional 
action is 

-L / d*+«V=s (r-2A- - J*) , (68) 



16tt 

where F q is the g-form field for stabilizing the g-sphere. The metric and the g-form flux are given by 

ds 2 = -di 2 + e 2/l *d£ 2 _i + P 2 dn 2 q , (69) 

and 

F q = be^...^ (70) 

where e^...^ is the volume element of the g-sphere. The Einstein and Maxwell equations reduce to two algebraic 
equations: 

(q-l)p- 2 -(p-l)h 2 =b 2 , 
{q-\) 2 p- 2 + {p-\) 2 h 2 = 2K. [ ' 



The entropy S is given by 



where $ is the total flux defined as $ = —bp q tt q and O p _2, g are respectively the volume of the unit (p— 2)-sphere and 
g-sphere. From (|71[) we can express S as a function of one parameter. In this case we can also see that the entropy 
S(<f>) as a function of <& is not single-valued and splits into a high-entropy branch and a low-entropy branch. The 
critical point dividing into two branches becomes 

a 2A(p-2) 2 _ (p-i)(p + q -2) a 

CTi (p-l)3(p + g _ 2 )' cri (p_ 2)( g -l) cri ' {<6) 

which is determined by dS/dh = 0. These values are nothing but the threshold at which a tachyonic mode appears 
in the scalar sector Q. Therefore, we have shown the equivalence of thermodynamic stability and dynamical stability 
in the general Freund-Rubin flux compactifications with dS p x S q . It is worthwhile investigating a similar relation in 
models with branes of codimensions higher than two. 
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APPENDIX A: THERMODYNAMIC RELATIONS 



In this appendix, we show the thermodynamic relations in the 6D braneworld solution with the 4D de Sitter space. 
The brane positions are determined by the two positive roots of f(r) = 0, where 



f(r) = ^r 2 - ^ - —. (Al) 
J\ > 10 r 3 12r 6 

Eliminating Ai from the conditions /(l) = and /(a) = 0, we obtain 



h 2 (l * 3 ) - ^(1 - « 5 ) - ^(a 3 - 1) = 0, (A2) 



and this yields 



k + + = + l/'(a)a 4 = ^(1 - a 5 ) - ^(a 3 - 1) 

= /i 2 (l-a 3 )--^(a 3 -l) (A3) 
= —.4-6$. 

47T 

This is the relation between the total area of the horizon A, the total magnetic flux of U(l) field <I>, and the proper 
periods at the brane positions K± which are respectively defined as 

$ = J F rip dr A d(f> = ^{a 3 - 1)A0, 

and 

K+ = -i/'(l)A& «_ = i/'(a)A^. (A5) 

Here A0 is a given period of the angular coordinate of extra dimensions. Recalling Eq. (j4|), we find k± = 2tt — &ita± 
where a± are tensions of the r± -brane. We note that Eq. (jA3|) is similar to the thermodynamic relation for the 6D 
RNdS black hole. 

Furthermore we can obtain the differential relation for the above quantities which is similar to the first law of black 
hole mechanics. The parameters of this spacetime are determined by two equations, /(l) = and f(a) = 0; then the 
variation of parameters in the two equations becomes 

f'(a)da + 2hdh - ^ - -^db = 0. 

a 6a (A6) 

2hdh -dM- -db = 0. 
6 



Eliminating dM, we have 



/,-. <ln ' - ,4<l ( ^- ) + <!><!/, = (). (AT) 



Although the period of the angular coordinate, A0, is necessary for determining the global geometry of the space- 
time, locally it is irrelevant. In other words, A, $ and k± in Eq. (|A3|) are "extensive" quantities with respect to the 
normalization of Acf>. Therefore we express the above relations in terms of following "intensive" quantities, 77 = k + /k_, 
A = A/k—, and $ = <5>/k_. Using these quantities normalized by k_ it is convenient to discuss the case with the 
tension of the r_-brane fixed. From Eqs. (|A3[) and (|A7|) . we have 
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and 



da 4 -^df — ] l , I>(l/ , : : „. 

47T 



Combining these two equations we obtain another expression, 

3h A - ~ 
dri = dL4 - 6d$. 

47T 



(A9) 



(A10) 



When the tension of the r_-brane is fixed, the variation of r\ is directly related to the variation of the tension of the 
r + -brane, i.e., dr\ = — drr + . In this sense the above expression is important. 

APPENDIX B: EUCLIDEAN ACTION 

In this appendix, we revisit the derivation of the thermodynamic relations from the Euclidean action. (See, for 
example, [ioj.) We take an ansatz for the Euclidean 6D metric with topology S 2 x S* 4 , 



ds 2 

and for the 17(1) field, 



X{r) 2 d<t> 2 + Y(r) 2 dr 2 + r 2 [(l - h 2 p 2 )dr 2 + (1 - tfp^dp 2 + p 2 dfl 2 ], 



A M dx M = A(r)dcj>, 



(Bl) 



(B2) 



The 4-dimensional part of the metric is a 4-sphere with the radius h, which is the Euclidean continuation of the 4D 
de Sitter space with the Hubble parameter h by the Wick rotation of the time coordinate r — > it. The Euclidean 
time t has a period At — 2n/h related with the de Sitter temperature. 
We obtain a component of the Einstein tensor, 



Y 2 



and the non-vanishing components of the energy-momentum tensor are 



=-(— V 

2 J 



The nontrivial Maxwell equation is 



which is integrated to yield 



XY 



XY 



= 0, 



-b. 



(B3) 



(B4) 



(B5) 



(B6) 



where b is the integration constant. 

The component of Einstein equation, G^a, — a + A 6 = is 



Y 2 



and it is easily integrated to obtain 



1 

Y 2 



h 2 - ^r 2 
10 



M 



12r 6 



(B7) 



(B8) 



where M. is the integration constant. The Ricci scalar and the field strength are 
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and 



1 p t?M N 



( — 

\XY 



(BIO) 



The Euclidean action is given by 



= -^~J d "^V9 (R 2A 6 - l -F MN F MN ^ 



16 

167T 



/■i/n /-i 

/ P 2 dW 

JO J a 



dr { -4XY 



1 

y2 



r 4 AT' 



2A 6 r 4 XY" 



XY 



6/i 3 

-A(6$ + K+ + a 4 K _), 



(Bll) 



where 



and 



A</> 



Y 



= K+, A0 



X' 



= — K_ 



$ EE A0(A| r= l - A| r=Q ). 



(B12) 



(B13) 



Note that k± and $ are the boundary data held fixed in the current action. Here, j3 is the inverse of the temperature 
of the de Sitter horizon, i.e., [3 = At. 

Two 3-branes are located at r = a, 1, so we require 



X(a) = X(l) = 0. 

The regularities at r = a, 1 lead to 1/F 2 | r=a ,i = 0, and we obtain 



.2 A 6 2 



12a 6 



^ ^ A 6 b 2 

= h - — - — - M = 0. 
a 3 10 12 



From these two conditions we can determine the two integration constants 6, M : 



M(a,h) = (l + a 3 )h 2 



A 6 (l-a s 



10(1 -« 3 ; 



-, b 2 (a,h) 



6A 6 a 3 (l-a 5 ) 
5(1 -a 3 ) 



I2a 3 h 2 , 



The reduced action /* for fixed n± and $ is 



I* ($, k+, k_; a, h) = -^-j (6$ + k+ + a 4 K_). 



(B14) 



(B15) 



(B16) 



(B17) 



Extremizing this action with respect to a and /i gives us 



or it f Ob , , \ „ 



(B18) 



96 
9a 



These two conditions are reduced to 

$/k_ = 4a 3 
where 5* is defined as 



^h 4 

k+ + K_a 4 = —S* - 6$, 

7T 



7r$ 96 7ra 3 <I> 
r^ 3 ^ = bhF' 



(B19) 



(B20) 
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We denote (a, h) satisfying the conditions (|B19|) for given k± and $ as a c ($, k+, k_) and h c (Q, k+, and we have 
a classical action 



J[$, K+, K_] = /*($, « + , k_; a c , fte) = - g^J (^c* + + a*K-) 

= a e , h c ) = «+, k_], 

where b e = b(a e , h ). Moreover, the differential relation is given by 

dS 



-dl= ^j(6 e d$ + d/t H 



a a dK-) 



or by using "intensive" variables S, $ and 77 we obtain the alternative expression 

- - 3h 4 - 



dS = -^(bd$ + dr)), 77 + a 4 



-S-b$. 



(B21) 



(B22) 



(B23) 



From the conditions (|B19j) we obtain the stationary points of I* , in other words, a series of classical equilibrium 
solutions. To see the (thermodynamic) stability of the solutions, we have to know which of the equilibrium solutions 
has the minimum of / (or the maximum of S) for given <I> and k±. For example, requiring one of the conditions, 
dl* /da = 0, we regard /* as an one-parameter function of h. (Note that $ and k± are fixed.) Then, we have 



dJ* 


dr 


dl* (da 


~dh 




da \dh 


d 2 I* 


d 2 r 


d 2 r ( 


dh 2 


dh 2 


dhda V 



4>,K_ 

da N 



dr 

'dh 1 



(B24) 



At a stationary point they become 

d/* 
dh 



= 0, 



d 2 r 



dh 2 



a c ,h c 



a c ,h c 



dS* 
~dh 



dS* (da 
da \ dh 



<t\K_ 



1 f^L 

h \~dh 



(B25) 



Hence, the condition for taking a minimum value of / (a maximum value of S) for fixed $ and k± is given by 

< 0. (B26) 



as 

dh 



APPENDIX C: GAUGE CHOICE 



1. Harmonics on dS4 



We give definitions of scalar, vector, and tensor harmonics on the 4D de Sitter space, g^ is the metric of the 4D de 
Sitter space with the Hubble parameter h, namely, R^ va p = h 2 [g^ a g v p — g^f}g va \, and V M is the covariant derivative 
associated with t^,„. 

The scalar harmonics Y satisfies 

V 2 F = m 2 Y. (CI) 

The vector harmonics V(t)^ satisfies 

[V 2 - 3h 2 }V (T ^ = m 2 V (I> , V"V (T )„ - 0, (C2) 

and V(l) m is defined as 



V(L), = d„Y. 



(C3) 
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The tensor harmonics h^)^ satisfies 

[V 2 - 2h 2 }T (T) , u = m 2 T {T) ^, VT (1> „ = T (T) % = 0, (C4) 

and T^xLjLLjYJjul/ is defined by 



m 2 



(C5) 



= 2v M v„y - — 5„ v y, 

1. Gauge choice 

We expand the perturbed metric in harmonics of the 4D dc Sitter space with the Hubble parameter h as 

Sg M Ndx M dx N =h ab Ydx a dx b 

+ 2(h (T)a V (T)ti + h (L)a V { ^)dx a dx» (C6) 
+ (^(T) ^(T)^ + ft(TL)T(TL)f*v + ft(LL)T(LL)/ui/ + ^(Y^Y^Ma^da; 1 ', 

where y, V(t,l)> an d T(t,tl.ll.y) are scalar, vector, and tensor harmonics, respectively. Greek indices /i, v, ■ ■ ■ run 
over the 4D de Sitter space and Latin indices a, b, ■ ■ ■ run over the 2D extra dimensions (r, (f). The coefficients /i a b, 
h(T)a> ^(L)a! ^(t)i /i(TL) > ^(ll) > an d ^(y) depend only on r assuming that the perturbations are axisymmetric. The 
perturbations of the U(l) gauge field can be expanded as 

5A M dx M = a a Ydx a + {a {T) V {T)ft + a (L) V^dx" , (C7) 

The infinitesimal coordinate transformation and U(l) gauge transformation are 

SgniN — > SgMN - C-zgMN, SA M ^ 5A M + Omx- ^Am- (C8) 

Also, the gauge parameters can be expanded as 

£ M dx M = (£( T ) + £(L) V^ds" + £ a Ydx a , X = X(Y)Y. (C9) 

The tensor- type component is gauge invariant: 

^(T) — * ^(T) = ^(T)- (CIO) 

The vector-type components transform as 

2 



fr(TL) — » ft(TL) = ft-(TL) ~ £(T), ^(T)r — » ^(T)r = ft(T)r _ £(T) + ~£(T)> 
^(T)0 — * ^(T)0 = ^(T)</>, «(T) — * S(T) = a (T), 



(Cll) 



where the prime denotes the derivative with respect to r. Choosing £(t) = ^(tl)j we have /1(tl) = 0- With this gauge 
choice we use new variables. Using a part of perturbed Einstein equation we obtain algebraically 

V)r = ^ (C12) 

where C is an arbitrary constant for m 2 = — 6/i 2 , or C = for m 2 ^ —6h 2 . When m 2 = — Qh 2 , this constant 
corresponds to the pure gauge mode for tensor perturbations because V M T( TL ) M ^ = [V 2 + 3/i 2 ]V( T )„ = and T^l)^ = 
2V m V(t) ai = 0. Hence we can choose C = 0. 
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The scalar-type components transform as 



(LL) - 


■* fr(LL) = fe(LL) " f(L), 


/l (Y ) 


(L)r ' 


, 2 

-> fy L)r = ft, (L)r - C(L) + ~£(L) 






f 

-y h — h — if — —f 




htf,tf, - 


-+ ft.^^ = h^j, — ff£ r , 


0(L) 


a r - 


^a r =a r + x{ Y) - A , 


a<f, 



m 



f 

h r< p — h r< f, — ^ + "T^0' 



a (L) = «(l) + X(Y) - ^4y, 



(C13) 



Choosing 



£(L) = fylX), 



h r <t,{p) 
f(p) ' 



£r = «(L)r - k(LL) + ~^(LL) j 



X(Y) 



-«(L)+^, 



(C14) 



where C" is an arbitrary constant, we have ^I(ll) — ^(L)r = ^r</> = O(l) = 0. With this gauge choice we redefine each 
component as follows: 



V) = 



-(*i+3$ 2 )/, 



l (L)0 = h (L)<t>, h rr = ($l + $ 2 )/f, 

a r = a r , »0 = a^. 



(C15) 



Moreover, using parts of the perturbed Einstein and Maxwell equations we algebraically obtain the following relations 
for four variables: 



* = $2, h (L)<j> = Cf, a r 



A' 1 

-y^(L)0 = -CA', = — 



/'* 2 + ^(/r 2 $i)' 
2r z 



(C16) 



where C is an arbitrary constant for m 2 — 0, or C — for m 2 ^ 0. This constant corresponds to the residual 
(homogeneous) gauge freedom, £^ = Cf and X(y) = CA. We can set C = 0, i.e., fy^u = a r = 0. 
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